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ABSTRACT 

We present a single-scattering formalism for incoherent resonant light scattering by dilute 
quantum gas systems such as the atomic-trap Bose-Einstein condensates. We show that 
resonant scattering gives access to more information than the dynamical structure factor, 
familiar from non-resonant scattering. In particular, we show that the detuning dependence 
of the incoherent scattering cross-section allows the direct determination of the BEC pairing 
density (ipip) , which is a broken symmetry and provides evidence that the condensate is not 
in a good number state. 
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The unusual properties of the atomic-trap Bose-Einstein condensates PJ— make them 
prime examples of dilute many-body systems with highly interesting microscopic structures. 
Experimentally, the most convenient probe is resonant optical scattering. From the theo- 
retical perspective, this poses an interesting problem - resonant scattering is a second-order 



as Javaneinen pointed out || ||, a Van-Hove-like expression is recovered in the off- resonant 
limit. In this paper, we discuss a single scattering formalism that is valid for arbitrary values 
of the detuning. 

Typically, the width 7 of the resonant transition is much larger than the trap frequency ujt 
(twice the ground state energy of the trapping potential, K = 1 in our units). Nevertheless, 
for some long-lived states, 7 can be comparable to the relevant excitation energies of the 
many-atom system, so that the scattering is fast on the scale of the excited atom motion, 
but not necessarily on the scale of the many-body dynamics. In that case, we find that 
resonant scattering gives access to more information about the many-body structure than 
the dynamical structure factor. In particular, for a dilute Bose condensate we are lead to the 
remarkable conclusion that resonant scattering allows a direct determination of the pairing 
density (ipip), related to the Bose symmetry breaking. 

It is instructive to recall the non-resonant scattering result. We adopt the convention 
that for an incident particle of momentum k; n and frequency u>i n , scattered into a state of 
momentum k out and frequency u) ont , q is the momentum transferred to the target system, 
q = k in — ko Ut , and uj = u>i n — u ont . Van Hove showed that the differential cross section 
d 2 a/dQ duj, where dQ is an infinitesimal solid angle, only depends on q and uj, and is equal 
to 



where /(q) is the scattering length describing the scattering of an incident particle by an 
individual target particle, and 5(q, uj) is the dynamical structure factor of the many-body 
system. The structure factor is the Fourier-transform of the density- density correlation 



process and cannot be described by means of the usual Van Hove theory [|J . Nevertheless, 
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function, 

S(q) = (2-k)- 1 J d 3 x d 4 x' exp[iq ■ (x - x')] (p(x')p(x)) , (2) 

where ( ) denotes the thermally averaged expectation value over the initial states of the 
many-body target system, and where we introduced the four- vector notation, q = (q, a;), 
x = (x, t), and q ■ x = q • x — out. In second quantization, the density operator is equal to : 

p(x) = ^(x)tJj(x) , (3) 

where ip and "ijy are the annihilation and creation fields in the Heisenberg picture. The 
integration / d 3 x in (||D is over the three spatial components, the time component being 
fixed. 

Unlike non-resonant scattering, resonant photon scattering involves two interactions - the 
atomic excitation while absorbing the incident photon, and the de-excitation accompanied by 
photon emission. In summing over scattering histories, we integrate over the photon absorp- 
tion amplitude at x\, the photon emission amplitude at X2, and the amplitude for the excited 
atom to move from X\ to X2, which is the excited atom propagator G(x2, where X\ 7^ x<i 
0. To deal with this type of non- locality, it is useful to work in the Wigner- representation, 
where a two-point function, /(xi,^) in coordinate representation, is represented by 

fw( x 'iP) — J d A r f(x + r/2, x — r/2) exp[ip ■ r] . (4) 

Similarly, we introduce the Wigner distribution operator: 

Pw{x;p) = (27r)~ 4 / d 4 r ffi (x + r/2)^(x — r/2) exp[ip • r] , (5) 



the expectation value of which gives the Wigner distribution [|J, 0, the quantum analogue 
of the classical phase-space distribution function. A straightforward but somewhat lengthy 
derivation yields the following cross-section for light scattering involving an atomic transition 
of natural frequency u>o and dipole moment d : 

d 2 a[e in ,e ont ;A] _ /3 7 \ 2 , v 2 1 r 

dQ dihuj) ~ \W (m )( ° ut } 2^ X J dxdX dpdp 
exp [iq ■ (x - x')] G w (x;p + k+) G* w (x'\p' + k+) (pw(x' ,p')p w (x,p)) , (6) 



where A is the detuning of the incident photons, A = oj\ n — u , k + is the average of the incident 
and outgoing photon momentum four-vector, k + = (k + , u + ) = ([k in + k ou t] /2, [u)in + uj out ] /2), 
d = d/|d| and e in and e out are the polarization directions of the incident and detected pho- 
tons. 

The propagator G can be computed from the eigenstates of the effective potential V'(x) 
experienced by the excited atoms. Nevertheless, if the excited state lifetime is so short 
that the excited atom experiences a change in V small compared to the width 7, (i.e. 

[F • v] /7 2 I < 1, where F is the force and v the velocity of the excited atom), then Gw is 
approximated accurately by the homogeneous propagator with V(x) as position dependent 
shift: 

G w (x;p + k + ) = ■ 1 — r 

p + uj + - [H e (p + k + , x) + uj \ 

= I (7) 

where po is the frequency-component of the momentum four-vector p, and H e is the self- 
energy of the excited atom, including the potential V and the width —27/2. Expanding the 
propagator Gw in powers of the inverse detuning, gives a A _1 -expansion of the cross-section, 
with coefficients that are n-th order moments of the [H e (p + k + ; x) — (po — u;/2)] -functions 
of the propagators (0) with respect to the (p(x'; p')p(x; p)) -correlation function. At first 
sight, the appearance of k + , instead of the momentum of the excited atom, might seem 
puzzling. The problem is resolved when we return to the ordinary coordinate representation 
using 



d^x d^p exp[iq ■ x] F(p) p w (x;p) = j d 4 x exp[iq ■ x] ^)\x) JF (p + q/2) t[j(x) 

d 4 x exp[iq ■ x] $(x) T (p — q/2) ip(x) 



where F(p) is an arbitrary function, and T{p) the operator obtained by replacing the p- 
vector in the expression for F, by the p-operator, pj = -rd/dxj, if j = l,2 or 3, andpo = id/dt. 
The right — ^ and left ^— arrows indicate that the p-operator only acts upon the field operator 
immediately to the right or left. Since q/2 + k + = k in , the [H e (p + k + ; x) — (p — uj/2)]- 
functions give an operator 7Y c (k in + p; x) — id/dt. The time derivative id/dt, acts on a 'hole' 



state. In this sense, 7i c — id/dt describes the evolution of the excited atom-hole pair and we 
denote the operator H e — id/dt by TC e - g - The cross-section is then equal to 



d 2 a [6 in , e ou t; A] 
dQ d(huo) 
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(37/4/c)(e in • d)(e out • d)* — x / d 3 x d 4 x' exp [iq ■ (x — x') 

2,71 
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i/j(x')if}' (x 



A- Hc- g (x 



7p{x)) 



(9) 



A— 7~t e —g 

For large detunings, we approximate the denominators in (|9|) by A, giving the off- resonant 
limit reported by Javaneinen [[J : 

2 

S(q,u), (10) 



d 2 a [e in ,e out ; A] 
A^oo dQ d(huj) 



(3 7 /4A;)(ei n -d)(e out -d) 



A 



where we recognize the off-resonant limit of the scattering length for optical resonant scat- 
tering from atoms, showing the similarity of (|i~0"D to the Van- Hove expression ([!]). 

The general result @ indicates how and when the off-resonant limit (|T0|) breaks down. 
Especially for cold fermionic atom systems, it can be important to correctly include the 
effects of the Doppler shifts, contained in (P). For the atomic Bose condensates, the recoil 
energies can be of importance, as we show below. 

A straightforward generalization of the customary detailed balance argument [|10| for the 
structure factor shows that changing the sign of the energy transfer to — > — uj (while keeping 
A constant) similarly gives : 



d 2 a , . d 2 a . . 

-(A; uj) = exp(puj) (A — uj; —uo) 



(11) 



dQ duo dQ duo 

where (3 is the inverse temperature. Verifying (|TTD experimentally (which does note re- 
quire angular resolution) for the atomic-trap BEC systems can test the hypothesis that the 
observed condensate is in thermal equilibrium and provides a direct measurement of the 
temperature. 

We now discuss scattering from a dilute homogeneous BEC at zero temperature. Since 
our concern is not with coherent scattering (observed at uo = , q = 0), we calculate the 
incoherent part of the cross-section, d 2 a nc /dQduo, subtracting (fo ip) (ip^ tp) from (ift'ipip^ip) in 
(|). We expand the -^-fields in plane wave states, tp(x) = V^ 1 ^ 2 a^it) exp [ik ■ x], fo(x) = 



V~ x l 2 Q>k(t) exp [— ik • x], where V is the volume (V — > oo in the end), and make the 
mean-field Bogoliubov approximation, treating ak=o, a k=0 as c-numbers, ak=o, a k=0 ^ V^ch 
keeping terms up to order N Q . The result is 

(far, 



dfl duj 



(3 7 /4A;)(e in • d) (e out ■ d) 5(u - c^ q ) x 



1 



+ 



+ 



A +h/2 - ujc, 
1 

A -27/2 - 
1 

A +27/2 - CJ q 



(aL q a_ q ) + 



1 



1 



A +ij/2_ 
1 

A -17/2 



A +«7/2 

(0-q°q) 



(a q o_ q ) 



(a q a q ) 
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where A is the effective detuning, A= A — £>(|k|) — V", and £>(|k|) represents the kinetic 
energy of an excited atom with momentum equal to the resonant wave number. In addition 
to the peak at u = co> q represented by fll2D , the incoherent scattering spectrum at fixed 
scattering angle and detuning, has a continuous background, caused by k-modes (k 7^ 0) 
excited to k + q states. In most cases of interest, the ratio of the integrated background 
intensity to the peak intensity is of the order of the depletion (N — N ) /N , where N is the 
total number of atoms. 

The Bogoliubov approximation implies that in each process a particle is taken out of, or 
put into the condensate. This leaves only two possibilities if q is the momentum transfer: 
a particle is excited from the condensate and ends up into a final state of momentum +q 
(corresponding to the subscript q in (|T2"D), or a particle is taken from the initial state of 
momentum — q and ends up in the condensate (corresponding to the — q-subscripts). In the 
+q-process, the excited atom has momentum k in , in the — q process, the excited atom has 
momentum — q + k; n = k ou t- I n either case, the excited atom has a center-of-mass energy 
that is u>( |k| ) + V. On the other hand, the initial-state energy for the q and — q processes 
are different, so that the 7i e - g -operator gives cD(|k|) + V — 27/2 for the +q-process and 
a>(|k|) + V — ujq — 27/2 for the — q scattering. The — q process leaves the condensate with 
an extra particle, whereas in the +q process, a particle is removed from the condensate. In 
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the usual single particle picture, one could expect these processes to give orthogonal final 
states, thereby precluding any interference. Nonetheless, the condensate is expected to be in 
a coherent state, rather than a number state, and the +q and — q final states have a finite 
overlap proportional to (a q a_ q ). Conversely, detecting the interference of the +q and — q 
scattering events answers the question wether or not the condensate is in a good number 
state. With regards to this isssue, which has received considerable attention in the recent 



literature |11|-|15|, we note that the incoherent scattering scheme provides an alternative to 
the interfering condensate experiments. 

To experimentally detect the interference, we note that the difference in recoil energy for 
the +q and — q scattering processes gives a distinct detuning dependence to the interference 
contribution in the cross section (|T^): 

d 2 a n 



dQ duo 
( 



tin; t ou t ; A] « iVo (37/4fc)(e in • d) (e out ■ d) 8{w — ^ q ) x 



^(A-Wq) +( 7 /2) 



(al q a_ q ) + (a q a_ q ) 



+ 



A) +(7/2) 2 
\ 



1 + (a q a q ) + (a q a_ q ) 



X 



A -u t 



(7/2) 2 



A 



(7/2) 5 



■(a q a_ q ) 



(13) 
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where we used that (aL q a_ q ) = (a q a q ), (a q a q ) = 1 + (a q a q ) and (a^L q a q ) = (a q a_ q ). 
Thus, the detuning dependence of the peak intensity differs from the simple Lorentzian 



1 -1 



A +( 7 /2) ; 



. The simplicity of the actual detuning dependence, which contains only 



two parameters, the occupation number (a q a q ) and pairing matrix element (a_ q a q ), sug- 
gests a simple fitting of the experimental curves to determine their values. In Fig.(l) we 
show the detuning dependence of the intensity in the peak of the incoherent scattering spec- 
trum, for the special case that the scattering angle corresponds to a momentum transfer for 
which the excitation energy c<j q is equal to the chemical potential fj, and \x is equal to the 
width 7. The full curve shows the actual intensity, whereas the dotted lines indicate the 
contributions proportional to the occupation number and the pairing matrix element. The 
occupation number and pairing matrix elements where calculated in the Bogoliubov approx- 
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imation at T=0 (the finite temperature generalization is straightforward, see for example 
fTE|). Requiring the non-Lorentzian contribution to be measurable leads to the condition 
that u q is of the order of 7. The Bogoliubov theory at T=0 gives (a q a_ q ) = — /i/2u; q , yield- 
ing a magnitude of the non-Lorentzian term relative to the other contributions that is of 
the order of ~ (t^X^)- Therefore, the best signal is obtained for backscattering, q = 2k, 
cj q ~ 4u r , where u r is the recoil energy, with a relative magnitude of ~ 2 (^) (^2)- 

Of course the atomic-trap system is not homogeneous. In fact, we cannot require the 
condensate to be too close to homogeneity because if the system is optically thick, the single 
scattering approximation breaks down. Nevertheless, it is possible to have an optically thin 
system for which fi ^> ut, in which case the Thomas- Fermi description is valid and the 
condensate |fT7| , [IS|, as well as the fluctuations JTP[, |SIJ behave locally in the same manner 
as the homogeneous system. In a Thomas-Fermi description of the incoherent scattering, the 



single peak is 'broadened', giving a feature in the frequency interval from y {q 2 /2m + /i) 2 — y? 
to q 2 /2m. Within this feature, the intensity of a frequency interval u' to u' + du has 
information about the spatial region in which the local excitation with momentum q has an 
excitation energy cj q (r) in the du - interval, u' < cj q (r) < u' + du. If the density is nearly 
constant in this spatial region, the detuning dependence of the intensity in the du - interval 
is described by the above theory. 

In conclusion, we have shown that resonant light scattering can measure the pairing 
matrix elements of Bose-Einstein condensates. The finite value of the 'pairing' density, (ipip), 
is an example of a higher-order broken symmetry and is of fundamental interest. Measuring 
its value would constitute a detailed test of mean-field theories and finding a non-zero result 
would prove experimentally that the condensate is not in a number state. 
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Figure Caption 



Figure 1: Plot of the detuning dependence of the peak intensity in the incoherent scat- 
tering spectrum. The calculation was performed in the Bogoliubov approximation at zero 
temperature for the special case that the momentum transfer corresponds to an excitation 
energy u; q equal to the chemical potential //, and /i is equal to the excited state width 7. 
The full curve shows the actual intensity, the dotted line with negative values shows the con- 
tribution proportional to the pairing number (a q a_ q ) (which is negative). The dotted line 
with the positive values represents the contribution proportional to the occupation number 
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